Entanglement renormalization 
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In the context of real-space renormalization group methods, we propose a novel scheme for quan- 
tum systems defined on a D-dimensional lattice. It is based on a coarse-graining transformation 
that attempts to reduce the amount of entanglement of a block of lattice sites before truncating its 
Hilbert space. Numerical simulations involving the ground state of a ID system at criticality show 
that the resulting coarse-grained site requires a Hilbert space dimension that does not grow with 
successive rescaling transformations. As a result we can address, in a quasi-exact way, tens of thou- 
sands of quantum spins with a computational effort that scales logarithmically in the system's size. 
The calculations unveil that ground state entanglement in extended quantum systems is organized 
in layers corresponding to different length scales. At a quantum critical point, each rellevant length 
scale makes an equivalent contribution to the entanglement of a block with the rest of the system. 

PACS numbers: 



Renormalization, one of the conceptual pillars of sta- 
tistical mechanics and quantum field theory, revolves 
around the idea of rescaling transformations of an ex- 
tended system (lj . These so-called renormalization group 
(RG) transformations are not only a key theoretical el- 
ement in the modern formulation of critical phenomena 
and phase transitions, but also the basis of important 
computational methods for many-body problems. 

In the case of quantum systems defined on a lattice, 
Wilson's real-space RG methods [2], based on the trun- 
cation of the Hilbert space, replaced Kadanoff's spin- 
blocking ideas [3j with a concise mathematical formu- 
lation and an explicit prescription to implement rescal- 
ing transformations. But it was not until the advent of 
White's density matrix renormalization group (DMRG) 
algorithm [J] that RG methods became the undisputed 
numerical approach for systems on a ID lattice. More re- 
cently such techniques have gained renewed momentum 
under the influence of quantum information science. By 
paying due attention to entanglement, algorithms to sim- 
ulate time-evolution in ID systems [B| and to address 2D 
systems @ have been put forward. 

The practical value of DMRG and related methods is 
unquestionable. And yet, they notably fail to satisfy 
one of the most natural requirements for a RG transfor- 
mation, namely to have scale invariant systems as fixed 
points. Instead, for such systems the number of degrees 
of freedom (that is, the Hilbert space dimension) of an ef- 
fective site increases with each rescaling transformation. 
This fact conflicts with the very spirit of renormalization, 
but it is not only disturbing from a conceptual viewpoint. 
The cost of a computation depends on the size of the 
effective sites and it becomes unaffordable after a suffi- 
ciently large number of iterations. 

In this Letter we propose a real-space RG transfor- 
mation for quantum systems on a D-dimensional lattice 
that, by renormalizing the amount of entanglement in 



the system, aims to eliminate the growth of the site's 
Hilbert space dimension along successive rescaling trans- 
formations. In particular, when applied to a scale invari- 
ant system the transformation is expected to produce a 
coarse-grained system identical to the original one. Nu- 
merical tests for critical systems in D = 1 spatial dimen- 
sions confirm this expectation. Our results also unveil 
the stratification of entanglement in extended quantum 
systems, and open the path to a very compact description 
of quantum criticality. 

Real-space RG and DMRG. As originally intro- 
duced by Wilson [2j , real-space RG methods truncate the 
local Hilbert space of a block of sites in order to reduce 
its degrees of freedom. Let us consider a system on a 
lattice C in D spatial dimensions and its Hilbert space 
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where s € C denotes the lattice sites and V s has finite 
dimension. Let us also consider a block B C C of neigh- 
boring sites, with corresponding Hilbert space 
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In an elementary RG transformation, the lattice £ is 
mapped into a new, effective lattice £ , where each site 
s' € £ is obtained from a block B of sites in C by coarse- 
graining. More specifically, the space Y' s , for site s' e £ 
corresponds to a subspace §g of Vg, 

Y' s , =§sC¥b, (3) 

as characterized by an isometric tensor w, Fig ([!]), 

w:Y' s ,^Y Bl vJw = I. (4) 



Isometry w can now be used to map a state |\&) € Yc 
[typically, a ground state] into a coarse-grained state 
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FIG. 1: Isometries and disent anglers. Left: a standard nu- 
merical RG transformation builds a coarse-grained site s', 
with Hilbert space dimension m, from a block of two sites si 
and S2 through the isometry w of Eq. ((4]). Right: by using the 
disentanglers iti and U2 of Eq. ([7]), short-range entanglement 
residing near the boundary of the block is eliminated before 
the coarse-graining step. As a result, the coarse-grained site 
s requires a smaller Hilbert space dimension m, m < m. 



I*') £ V&- A thoughtful selection of subspace §g is es- 
sential. On the one hand, its dimension m should be as 
small as possible, because the cost of subsequent tasks, 
such as computing expectation values for local observ- 
ables, grows polynomially with m. On the other, 
needs to be large enough that retains all relevant 
properties of White identified the optimal choice 

as part of his DMRG algorithm Q. Let p 1 ^ denote the 
reduced density matrix of |\&) on block B. Then, 
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(5) 



where \pi) are the to eigenvectors of p^ B ' with largest 
eigenvalues pi, and to is such that e > 1 — Xa=i Pi, where 
e is a preestablished truncation error, £-Cl. 

Notice that to refers to the amount of entanglement be- 
tween B and the rest of the lattice, C—B, as characterized 
by the rank of the truncated Schmidt decomposition 
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In other words, the performance of DMRG-based meth- 
ods depends on the amount of entanglement in 

Entanglement Renormalization. We propose a 
technique to reduce the amount of entanglement between 
the block B and the rest of the lattice C while still ob- 
taining a quasi-exact description of the state of the sys- 
tem. This is achieved by deforming, by means of a uni- 
tary transformation, the boundaries of the block B before 
truncating its Hilbert space, see Fig. ([I]). 

Let us specialize, for simplicity, to a ID lattice and to a 
block B made of just two contiguous sites s\ and S2. Let 
rq and r 2 be the two sites immediately to the left and to 
the right of B. Then we consider unitary transformations 
u\ and u 2l the disentanglers, acting on the pairs of sites 



nsi and s 2 r 2 , 

Ui : V n <g> " 
u 2 : V S2 (g) ' 



V ri (8 V Sl , u\u\ = u\u\ = I, 
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\u 2 = u 2 u\ = I. (7) 



Properly chosen disentanglers reduce the short-range en- 
tanglement between the block B and its immediate neigh- 
borhood. The original state of block B, 



.,[«] = 



(8) 



is now replaced with a partially disentangled state p^ , 

= tr riI . 2 [( Ul ® u 2 )p [rislS2r2] (m ® uajt] , (9) 

which has a smaller effective rank to, to < to. Our RG 
transformation consists of two steps: 

(?) First we decrease or renormalize the amount of en- 
tanglement between block B and the rest of C. This is 
achieved with disentanglers that act locally around the 
boundary of B. The block does not become completely 
disentangled because only entanglement localized near its 
boundaries can be removed. 

(ii) Then, as in Wilson's proposal, we truncate the 
Hilbert space of block B, and we do so following White's 
idea to target the support of the block's density matrix. 
But instead of keeping the support of the original p^ , we 
retain the (smaller) support of the partially disentangled 
density matrix p^ . 





FIG. 2: Left: MERA for a ID lattice with periodic boundary 
conditions. Notice the fractal nature of the tensor network. 
Translational symmetry and scale invariance can be naturally 
incorporated, substantially reducing the computational com- 
plexity of the numerical simulations. Right: building block 
of a MERA for a 2D lattice. Disentanglers and isometries 
address one of the x and y spatial directions at a time. 

The above steps, characterized by unitary and iso- 
metric tensors u and w, produce an alternative coarse- 
grained lattice C and a coarse-grained state |^) € 
that contains less entanglement than \ty') S V/v. Impor- 
tantly, the expectation value (\f r |o|\I') = tr [opf 7 ^] , where 
o is an observable defined on a small set of sites 1Z C C, 
can be efficiently computed from just \^f), u, w and o. 
There are two ways: (i) as described in [9j, p^ can be 
computed from fP 1 ^ = tr^_ 7 j|5')(^|, where TZ is the set 
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FIG. 3: Scaling of the entropy of entanglement in ID quan- 
tum Ising model with transverse magnetic field. Up: in a 
critical lattice [h — 1 in Eq. (|fOp ]. the unrenormalized entan- 
glement of the block scales with the block size L according to 
Eq. Instead, renormalized entanglement remains con- 

stant along successive RG transformations, as a clear man- 
ifestation of scale invariance. Line (i) corresponds to using 
disentanglers only in the first RG transformation. Line (ii) 
corresponds to using disentanglers only in the first and second 
RG transformation. Down: in a noncritical lattice [h = 1.001 
in Eq. (|10[> ]. the unrenormalized entanglement scales roughly 
as in the critical case until it saturates (a) for block sizes com- 
parable to the correlation length. Beyond that length scale, 
the renormalized entanglement vanishes (b) and the system 
becomes effectively unentangled. 



of sites of C causally connected to 1Z through u and w; 
(ii) alternatively, we can use u and w to lift o from C to 
£, producing o, and exploit that (#|o|\&) = (*|6|*>. 

The lifting of linear operators from £ to C generates, 
by iteration, a well-defined RG flow in the space of hamil- 
tonians such that, in the ID case, an interaction term h 
acting on at most three consecutive sites of C is mapped 
into an interaction term h acting also on at most three 
consecutive sites in C [analogous rules apply in D > 1 
spatial dimensions]. Thus, in spite of the fact that disen- 
tanglers deform the original tensor product structure of 
VjC, the above rescaling transformation preserves locality, 
in that it maps local theories into local theories. 

In the same way as DMRG is related to matrix product 
states (MPS) [7[, the above RG transformation is natu- 
rally associated to a new ansatz for quantum many-body 
states on a D dimensional lattice. The multi-scale en- 
tanglement renormalization ansatz (MERA) consists of a 
network of isometric tensors (namely the isometries u and 
disentanglers w corresponding to successive iterations of 
the RG transformation) locally connected in D + 1 di- 



mensions. The extra direction r, related to the RG flow, 
grows only as the logarithm of the lattice dimensions, see 
Fig. $Z§. Several properties of the MERA together with 
its connection to quantum circuits are described in Q . 
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FIG. 4: Spectrum of the reduced density matrix of a spin 
block. Up: as the size L of the spin increases, the number m 
of eigenvalues {pi} required to achieve a given accuracy e, see 
Eq. (J5J, also increases. In particular, m grows roughly ex- 
ponentially in the number r = log 2 L of RG transformations. 
The spectrum resulting from applying disentanglers leads to 
a significantly smaller rh invariant along successive RG trans- 
formations. Down: spectrum of the reduced density matrix of 
2 T spins immediately before and after using the disentanglers 
at the r th coarse-graining step. These spectra are essentially 
independent of the value of r = 1, ■ ■ • ,14. 

Example.— Figs. ©-{II study the ground state of the 
ID quantum Ising model with transverse magnetic field, 



H = 



E ■ 

<«1,«2> 



(10) 



for an infinite lattice [1 01 ] . The entanglement entropy be- 
tween a block B of L adjacent spins and the rest of the 
lattice, defined in terms of the eigenvalues {p.;} of pl 8 l as 
S(B) = — J2iPi l°g2^ij scales at criticality (h = 1) with 



the block size L as [11| 



Sl ~ r log 2 L, 



(11) 



On the other hand, off criticality S(B) grows monoton- 
ically with L until it reaches a saturation value for a 
block size comparable to the correlation length in the 
system. Fig. ((3|) shows the effect of disentanglers on 
the entanglement entropy of a block. Most notably, the 
renormalized entanglement of a block in the critical case 
is reduced to a small constant value along the RG flow. 
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Correspondingly, as Fig. (J4J) shows, the Hilbert space di- 
mension m of a block of size L = 2 T can be kept constant 
as we increase t, allowing in principle for an arbitrary 
number of iterations of the scale transformation. Our 
calculations, involving the reduced density matrix of up 
to L = 2 14 = 16, 384 spins, have been conducted with 
Hilbert spaces of dimension m = 8, while keeping the 
truncation error e at each step fixed below 5 x 10~ 7 . We 
estimate that without disentanglers an equivalent error e 
requires a Hilbert space of m k 500 — 1000 dimensions 
for the largest spin blocks. This exemplifies the compu- 
tational advantatges of using disentanglers. 

Discussion.— An appealing picture emerges from 
the above numerical exploration. Entanglement in the 
ground state of the quantum spin chain is organized in 
layers corresponding to different length scales in the sys- 
tem. The entanglement of a given length scale can be 
modified by means of a discntangler that acts on a re- 
gion with linear size as large as that length scale. We 
obserbe two situations: (i) Off criticality, the entangle- 
ment between a block of sites B and the rest of the lattice 
C consists roughly of contributions from layers of length 
scale not greater than the correlation length in the sys- 
tem. Therefore, after a sufficiently large number of RG 
transformations with disentanglers, the effective ground 
state becomes a product state, that is, completely dis- 
entangled, (ii) At criticality, instead, the entanglement 
between 1Z and the rest of C receives equivalent contri- 
butions from all length scales (smaller than the size of 
B) - see @] for a justification of the logarithmic scaling 
of Eq. (fTT|) . In this case, by applying disentanglers we 
obtain a coarse-grained lattice identical to the original 
one, including an effective Hamiltonian identical (up to 
a proportionality constant) to the original one [12j |. In 
either fixed point of the RG transformation is 

attained after sufficiently many iterations. Non-critical 
theories collapse into the trivial fixed point (a product 
state) while critical theories correspond to a non-trivial 
fixed point (an entangled state expected to depends on 
the universality class of critical theory). We conjecture 
that this picture holds also for quantum lattices in D > 1 
spatial dimensions. Interestingly, a third type of fixed 
point is possible, in the case of a ground state with ei- 
ther long-range order or topological order [l3j |. 

In summary, we have presented a quasi-exact real- 
space RG transformation that, when tested in ID sys- 
tems, produces effective sites of bounded dimension and 
have (scale invariant) critical systems as its non-trivial 
fixed points. The key feature of our approach is a lo- 
cal deformation of the tensor product structure of the 
Hilbert space of the lattice, that identifies and factorizes 
out those local degrees of freedom that are uncorrelated 
from the rest. This deformation is such that local opera- 
tors are mapped into local operators, so that relevant ex- 



pectation values for the original system can be efficiently 
computed from the coarse-grained system. 

We conclude with two remarks. First, a scale invari- 
ant critical ground state leads to disentanglers u and 
isometries w that are the same at each iteration of the 
RG transformation. A single pair (u,w), depending on 
0(rh 4 ) parameters, is thus seen to specify the ground 
state, leading to an extremely compact characterization 
that deserves further study. Second, renormalization 
group techniques are extensively used also in classical 
problems, where the MERA can efficiently represent par- 
tition functions of lattice systems, extending the scope of 
this work beyond the study quantum systems. 
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